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FOREWORD 


This  report  was  prepared  by  the  Plastics  and  Composites  Branch, 
Nonmetallic  Materials  Division,  Air  Force  Materials  Laboratory.  The 
work  was  initiated  under  Project  7340,  "Nonmetallic  and  Composite  4 

Materials."  Task  734003,  "Structural  Plastics  and  Composites, 
was  administered  under  the  direction  of  the  Air  Force  Materials 
Laboratory,  Air  Force  Systems  Command,  Wright- Patter  son  Air  For- 
Base,  Ohio;  Dr.  J.  M.  Whitney  of  the  Plastics  and  Composites  Branch 
(LNC)  was  the  Project  Engineer. 

The  report  covers  work  conducted  January  1971  through  March  1972. 

This  technical  report  has  been  reviewed  and  is  approved. 
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T.  y .  REINHART,  Jfc. 

Acting  Chief,  Plastics  and 
Composites  Branch 

Nonmetallic  Materials  Division 
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ABSTRACT 


Extensive  use  of  energy  methods  in  conjunction  with  classical 
beam  mode  functions  has  been  used  to  obtain  approximate  solutions  to 
homogeneous,  or  symmetrically  laminated,  anisotropic  plate  problems. 
Because  of  the  existence  of  cross -elasticity  bending  stiffness  terms, 
the  beam  functions  do  not  satisfy  the  natural  boundary  conditions.  As 
a  result,  bending  moments  and  stresses,  which  are  of  practical  interest, 
may  converge  to  the  wrong  solution  or  may  not  converge  at  all.  Further¬ 
more,  bending  deflections,  buckling  loads,  and  fundamental  vibration 
frequencies  converge  very  slowly  for  highly  anisotropic  materials.  This 
report  shows  that  improved  results  can  be  obtained  for  anisotropic  plates 
which  contain  strong  cross-elasticity  effects  by  using  a  classical  Fourier 
analysis  which  satisfies  both  the  geometric  and  natural  boundary 
conditions.  Numerical  results  are  presented  for  bending  under  transverse 
load,  buckling  under  biaxial  compression  and  pure  shear,  and  natural 
frequencies  of  flexural  vibration.  Both  homogeneous  and  laminated 
plates  are  discussed. 
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Plate  length  and  width 

Elements  of  anisotropic  bending  stiffness  matrix, 

Q. .h^/12  for  homogeneous  plates 
ij 

Plate  thickness 

Distance  from  plate  mid-plane  to  interface  between  the 
k  and  k  +  1  layers  of  a  laminate 

Number  of  layers  in  a  laminated  plate 

Distributed  bending  moments 

Distributed  twisting  moments 

Distributed  normal  stress  resultants 

Distributed  shear  stress  resultant 

Integral  of  the  density  through  the  plate  thickness 

Distributed  normal  load  over  plate  surface 

Elements  of  the  anisotropic  reduced  stiffness  matrix 

Reduced  stiffnesses  rleative  to  x^y1  axes 

Time 

Plate  deflection 
Circular  frequency 
Normal  stress 
Shear  stress 

Angle  of  fiber  axis,  x' .  relative  to  x  axis  of  plate 


IV 


INTRODUCTION 


With  the  development  of  high  performs.nce  fiber  reinforced  composite 
materials  for  structural  applications  has  come  an  increased  interest  in 
solutions  to  anisotropic  plate  problems.  A  large  number  of  solutions  exist 
for  bending,  buckling,  and  free  vibration  of  specially  orthotropic  rectangular 
plates  in  which  the  principal  elastic  axes  are  parallel  to  the  sides  of  the 
plate.  Many  of  these  solutions  are  summarized  in  References  [1]  -  [3], 

In  most  structural  applications,  however,  fiber  reinforced  composites  are 
constructed  of  unidirectional  plies  in  which  the  fiber  axis  is  oriented  at  an 
angle  0  to  the  x  axis  as  illustrated  in  Figure  1.  For  such  a  composite 
symmetrically  laminated  about  th'i  mid -plane,  the  bending  response  is 
governed  by  the  flexural  equation  of  a  homogeneous  anisotropic  plate  [4,5], 
including  the  cross-elasticity  bending  stiffness  terms  and  D^. 

Energy  methods  have  been  used  [  6,  7,  8,  9,  10]  in  conjunction  with 
classical  beam  mode  functions  to  obtain  approximate  solutions  for  bending, 
buckling,  and  vibration  of  homogeneous,  or  symmetrically  laminated, 
anisotropic  plates  with  various  boundary  conditions.  Extensive  numerical 
results  for  rectangular  plates  appear  in  Reference  [  3].  Because  of  the 
existence  of  cross-elasticity  bending  stiffness  terms,  however,  the  beam 
functions  do  not  satisfy  the  natural  boundary  conditions.  Recent  work  by 
Fraser  and  Miller  [11]  involved  the  use  of  a  generalized  Ritz  method  in 
conjunction  with  a  Fourier  series  and  Lagrange  multiplier  technique  to 
obtain  an  upper  and  lower  bound  on  buckling  loads  for  homogeneous  aniso¬ 
tropic  plates.  Again,  as  in  the  previously  cited  literature,  the  natural 
boundary  conditions  were  not  satisfied. 


A«hton  [  10]  has  shown  that  unless  the  natural  boundary  conditions  are 
satisfied,  bending  moments,  shear  resultants,  and  edge  reactions,  which 
are  of  particular  interest  in  bending  problems,  m?y  converge  to  the  wrong 
solution  or  may  not  converge  at  all.  Furthermore,  for  highly  anisotropic 
materials  bending  deflections,  buckling  loads,  and  free  vibration  frequencies 
converge  very  slowly. 

ft  is  the  purpose  of  the  present  paper  to  show  that  improved  results 
compared  to  existing  energy  solutions  can  be  obtained  for  anisotropic  plates 
having  strong  cross -elasticity  effects  by  using  a  classical  Fourier  analysis 
which  satisfies  both  the  geometric  and  natural  boundary  conditions.  An 
exact  solution  to  the  governing  equations  is  obtained  in  the  form  of  a 
double  sine  series  plus  a  double  cosii.e  seilos.  A  procedure  similar  to 
that  used  by  Green  [  12]  and  by  Fletcher  and  Thorne  [  13]  on  isotropic 
plates  is  employed  to  exactly  satisfy  both  the  geometric  and  natural  boundary 
conditions.  Although  the  Fourier  method  is  applicable  to  any  of  the  classical 
boundary  conditions,  the  present  work  considers  plates  which  are  either 
simply- supported  on  all  edges,  or  simply -supported  on  two  opposite  edges 
with  the  adjacent  edges  clamped.  These  boundary  conditions  are  sufficient 
to  show  improved  results  compared  to  existing  energy  solutions.  The 
Fourier  approach  has  been  recently  applied  to  anisotropic  rectangular  plates 
having  all  edges  clamped  [  14].  Since  all  of  the  proper  boundary  conditions 
were  geometric  in  nature,  excellent  agreement  with  existing  Ritz  solutions 
was  obtained.  This  procedure  has  also  beer,  applied  to  un symmetrically 
larrr'-^tcd  anisotropic  plates.  [  15,16].  The  laminating  sequence,  however, 
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was  such  that  the  and  cross-elasticity  bending  stiffness  terms 
vanished.  This  simplification  lead  to  displacement  and  stress  function 
solutions  whi:h  could  be  expressed  in  terms  of  a  single  function  in  the  form 
of  a  double  trigometric  series. 


ANALYSIS 

For  the  purpose  of  completeness  it  is  appropriate  to  briefly  review  the 
Fourier  method  of  solution  as  applied  to  the  anisotropic  plate  equation. 
Denoting  partial  differentiation  by  a  comma,  the  governing  equation  for  the 
bending  and  flexural  vibration  of  a  homogeneous,  or  symmetrically 
laminated,  plate  with  pre -buckling  inplane  loads  included  is  of  the 
form  [  1,  2,  3] 


D  w,  +  4D  ,  w,  +  2(D  +  2D//)v/» 

11  xxxx  16  xxxy  12  66  xyyy 


+  4D  ,  w,  +  D  w,  +  P  w 

26  xyyy  22  yyyy  tt 


(1) 


=  N  >y,  +  2N  w,  +  N  w, 

x  xx  xy  xy  y  yy 


For  laminates  symmetrically  oriented  about  their  mid-pJane  [  3,4,5] 
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ij 


A 
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Q..  z2dz 
ij 


(2) 


where  K  =  L/2  for  a  laminate  having  an  even  number  of  layers  and 
(L  +  1)  /2  for  a  laminate  having  an  odd  number  of  layers. 


Bending  Under  Transverse  Load 


Consider  a  rectangular  plate  which  lies  in  the  region  0  <_  x  <_a,  0  <^y  <b 
(Fig.  1)  and  subjected  to  a  general  transverse  load  q  which  can  be  presented 
by  the  Fourier  series 


V*  V"  .  imrx  .  n-rr  y 

(x,y)  =  Z  Z  vsinvsab 


m=  1  n=  1 


(0  <  x  <  a,  0  <  y  <  b) 


A  solution  to  eq.  (1)^  in  the  absence  of  inplane  forces  and  inertia,  of  the  form 


w(x,y)  =  Wj (x,  y)  +  w2(x,y) 


is  sought  in  the  region  0<x<a,  0<y<b  where  and  independently 
satisfy  the  geometric  boundary  conditions.  It  is  assumed  that 


00  00 


'.  =  2  I 


.  rmrx  .  mr 

A  sin -  sin— r- 

mn  a  b 


m=  1  n=  1 


■<■2  i 


_  m-rr  x  niT 

B  cos - cos— r- 

mn  a  b 


m=  0  n=  0 


Differentiating  Wj  term-by-term  with  respect  to  x  twice  yields 


W  W  £ 

V1  T— ^  m  TT 

Z  Z  a2 

m=l  n=l 


.  .  mux  .  nTr ' 

A  sm - sin  — 

mn  a  b 


(0  <  x  <  a,  0  <_y  <_b) 
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Since  eq.  (6)  is  not  valid  on  the  edges  x  =  0,  x  =  a,  further  differentiation 
with  respect  to  x  cannot  be  accomplished  term-by-term  [  12}.  Assuming 
W1  xxx  Can  be  rePresented  by  a  cosine-sine  series,  partial  integration  leads 
to  the  result 


CO 


+  0  b 
rm  n 


miT  x  .  niry 

cos -  sin  , 

a  b 


where 


n 


rb 

=  rt"  J  fWl.xx,a'y)  -  W-«|U'y)]  dy 


n 


4 

ab 


J  fWl,xx(a-y)  +  W>XX(0>^1  sin2bX  d>' 


a  =  1,  m  even 
m 

-  0,  m  odd 

0  =  0,  m  even 

m 

=  1 ,  m  odd 


A  similar  procedure  applied  to  w,  ,  w_  ,  w_  ,  w_  ,  and  w_ 

™  l,yy  2,x  2,  y  2,  xxx’  2,  yyy 

leads  to  ten  more  sets  of  constants:  c  ,  d  associated  with  w  along  the 

mm  1 , yy 

edges  y  =  0,  b;  e  ,  f  associated  with  w_  along  the  edges  x  =  0,  a;  g  ,  h 

n  n  2,x  m  ni 

associated  with  c.long  the  edges  y  =  0,  b;  i  ,  j  associated  with 

2,y  oy  n  Jn 
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w_  along  the  edges  x  =  0,  a;  and  k  ,  1  associated  with  w_  along 
2,  xxx  mm  2,  yyy 

y  =  0,  b.  All  other  desired  derivatives  can  be  obtained  through  term-by-term 

differentiation.  Substituting  the  appropriate  derivatives  of  and  w^  along 

with  eq.  (3)  into  eq.  (1)  ,  equating  like  Fourier  coefficients,  and  solving  the 

resulting  algebraic  equations  yields  A  ,  B  _,  B_  ,  and  B  in  terms  of 
~  mn  mu  On  mn 

the  coefficients  a,b,c  ,...,1  The  latter  set  of  coefficients  along 
n  n  m  m 

with  are  determined  such  that  the  boundary  conditions  are  satisfied. 
Equation  (1)  also  yields  a  relationship  between  i  and  k  . 

For  a  plate  having  all  edges  simply- sup  prated  the  following 
classical  conditions,  designated  SS,are  applicable, 
at  x  =  0  and  a: 


w  =  M  =  -  D  w,  -  2D  ,  w,  -  D  w,  =0 
x  1 1  XX  1 6  xy  12  yy 


(8) 


at  y  =  0  and  b: 


w=M  =  -  D  w,  -  2D  w,  -  D  w,  =0 
y  12  xx  26  xy  22  yy 


(V) 


The  moment,  or  natural,  boundary  conditions  lead  to  the  following 
relationships 


Dnb  ,  Dnb 

e  -  — -  a  ,  f 


b 


n  ~  2D^nTr  n'  n  2D^mr  n 


n  i  0 


(10) 


°22Rb  ,  D22Rb  . 
6m  2D  m,  m  m  ZD^rmr  m 


m  i  0 


(11) 


Thus  eqs.  (8)  and  (9)  yield  an  infinite  set  of  simultaneous  equations  for 

the  constants  Brtrt,  a  ,  b  ,  c  ,  d  ,  e  .  f  ,  g_,  h  i  ,  j  ,1:  ,  and  1  . 

00  n  n  m  m  0  0  60  0  n  Jn  m  m 


Truncation  of  the  system  at  .m  =  k*  and  n  =  N  leads  to  a  4(M  +  N  +  2)  system 


I 

I 

I 


of  equations.  This  system  can  be  separated  into  four  sets  of  equations  by 
ranging  nr.  and  n  over  even  or  odd  integers.  Any  desired  degree  of 
convergence  can  be  obtained  by  taking  appropriate  values  of  M  and  N. 

For  the  purposes  of  this  paper  it  is  useful  to  also  consider  the  combined 
boundary  conditions,  designated  CS,  where  the  edges  x  =  0,  a  are  clamped 
and  the  edges  y  =  0,  b  are  simply -supported.  Thus,  at  x  =  0  and  b: 


w  =  w,  =0 
x 


(12) 


and  the  conditions  along  the  edges  y  =  0  and  b  are  those  of  eq.  (9) 


The  slope  w  vanishes  along  the  edges  x  =  0  and  b  if, 
1 ,  x 
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a  mA 
m  mn 


=  0,  >  Z  mA  =  0 

/ .  m  mn 


(13) 


m 


m 


for  all  n  =  1 , 2,  .  .  .  .  , 


while  w  vanishes  along  the  same  edges  if 
cf  x 


e  =  f  =  0 
n  n 


(n  =  0,1,2,.  .  .,) 


(14) 


The  previous  system  of  eqs.  for  SS  conditions  can  now  be  applied  to  the 
CS  case  by  replacing  eqs.  (10)  and  the  moment  portion  of  eq.  (8)  with  (14) 
and  (13)  ,  respectively. 

Using  contracted  notation  the  plate  stresses  and  moment  resultants 
are  calculated  form  the  constitutive  relations 
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a  «  «  zQ.W  k. 
1  1J  1 


(M  =  1,2,6) 


M.  =  D..  k. 
i  ij  ) 


where  o',  =  <r  ,  <r  =  <r  ,  <r,  =  <r  and  the  curvatures,  k.,  and  moments, 
1  x  2  y  6  xy  i 

M.,  are  defined  in  an  analogous  manner.  The  curvature-displacement 
relations  are  those  of  classical  homogneous  plate  theory 


k  -  -  w,  ,  k  =  -  w,  ,  k  =  -  2w, 
x  XX  y  yy  xy  xy 


Buckling  and  Flexural  Vibration 


Consider  the  case  of  uniform  membrane  loading  N  ,  i.  e.  , 


N  =  -  K  N  ,  N  =  -  K  N  ,  N  =  K  N 
x  10  y  20  xy  30 


N  >0 


where  K. ,  K  and  K.  are  known  constants.  Solutions  to  eq.  (1)  with 

1  U  O 


q  =  0  are  of  the  form 


w  =  W(x,  y)  e 


where  co  is  the  natural  vibration  frequency  and  t  denotes  time.  Substituting 
eq.  ^18)  into  eq.  (1)  and  taking  eq.  (17)  into  account  yields 


D..W,  +  4D  , W,  +  2(D  „  +  2D,,)  W  +  4D..W, 

11  xxxx  16  xxxy  1  l  66  xxyy  26  xyyy 


+  D00W,  -  P(o  W  •  K  N  W  +  2K  N  W  +  K,N_W,  =  0 

22  yyyy  1  0  xx  3  0  xy  3  0  yy 
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Solutions  to  eq.  (19)  can  be  obtained  by  assuming 


W(x,y)  =  Wx  (x,y)  +  W2(x,y) 

where  W,  and  are  of  the  same  form  as  w^  and  w^,  respectively. 
Following  the  same  procedure  as  in  the  previous  section  leads  to  four  sets 
of  homogeneous  equations.  Each  set  corresponds  to  a  different  mode  shape. 
Natural  vibrat;^n  frequencies  are  obtained  by  letting  =  0,  while  static 
buckling  becomes  the  special  case  of  vanishing  co.  Proper  values  of 
and  u)  are  determined  by  allowing  the  determinant  of  the  coefficient  mat/ix 
in  each  of  the  four  groups  of  equations  to  vanish. 


COMPARISON  TO  EXACT  SOLUTION 

An  exact  solution  is  available  [17]  for  the  bending  deflection  of  a 
simply- supported  anisotropic  plate  subjected  to  a  uniformly  distributed 
load  and  having  the  following  elastic  stiffness  properties 


Q11  (Q12  +  2Q66*  _  Q 
Q--  ^  =  ^’Q 


16 


Q 


26 


22 


22 


Q 


=  -  0.5 


22 


(20) 


This  problem  has  been  solved  by  the  Fourier  series  method  (FS)  with 
M  =  N  =  1,3,5,7,9,11,  and  13.  Comparison  of  the  maximum  deflection 
for  increasing  terms  with  the  exact  solution  (ES)  is  shown  in  Fig.  2.  The 
maximum  deflection  is  available  from  Reference  [10]  for  the  Ritz  method 
(RS)  ,  and  is  also  shown  in  Fig.  2  for  M  =  N  -  1, 3,  5,  7,  and  M  =  9,  N  =  8 
Since  the  Fourier  series  method  satisfies  all  of  the  required  boundary 
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conditions,  more  rapid  convergence  is  obtained  compared  to  the  energy- 
approach  in  which  only  the  geometric  boundary  conditions  are  satisfied. 

It  should  be  noted,  however,  that  direct  comparison  of  the  Fourier  method 
and  the  Ritz  method  is  difficult  as  given  values  of  M  and  N  lead  to  a  different 
number  of  non-vanishing  terms  in  the  series  representation  of  the  deflection 
for  the  two  approaches.  Because  of  the  symmetry  of  a  uniform  load,  for 
example,  M  =  N  =  7  yields  25  non-zero  terms  in  the  Ritz  method,  while 
only  16  non-zero  terms  appear  in  the  Fourier  solution. 

DISCUSSION 

In  the  work  by  Ashton  [  10],  it  was  shown  that  the  rate  of  convergence 
for  simply -supported  anisotropic  plates  employing  the  Ritz  method  in 
conjunction  ^ith  a  double  sine  series  was  dependent  on  the  anisotropy  ratio 
Q*  / Q! '  Thus,  for  highly  anisotropic  materials  such  as  graphite -epoxy 

XX  £  C* 

the  energy  approach  will  provide  a  very  slow  converging  solution  for 
simple -support  boundary  conditions.  This  is  illustrated  in  Fig.  3  for  a 
single -ply  composite  subjected  to  a  uniform  transverse  load  and  having 
the  following  properties  with  respect  to  the  material  symmetry  axes 


These  are  typical  properties  of  a  high  modulus  graphite -epoxy  composite. 
The  maximum  deflection  is  shown  for  a  square  plate  as  a  function  of  fiber 
orientation  0.  Fourier  series  results  are  based  on  M  =  N  =  13,  while 
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the  Ritz  solution  is  based  on  M  =  N  =  7.  In  the  case  of  the  energy  method 

this  is  the  maximum  number  of  terms  available  in  the  work  of  Ashton  and 

Waddoups  [  6] .  The  difference  between  the  Fourier  solution  and  the  energy 

solution  is  a  maximum  at  orientations  for  which  the  Dw  and  D_.  cross- 

16  Z6 

elasticity  terms  have  the  greatest  effect.  Convergence  of  the  energy  method 

is  immediately  improved  when  two  sides  are  clamped  as  also  illustrated  in 

Fig.  3.  This  improvement  is  due  to  the  fact  that  the  boundary  conditions 

are  all  satisfied  on  the  clamped  side.  For  all  sides  rigidly  clamped  the 

Fourier  method  and  the  energy  method  both  converge  rapidly  [  14]. 

In  Fig.  4,  the  bending  moment  across  the  centerline  y  =  b/2  is  shown 

for  the  45°  oriented  simply- supported  plate  of  Fig.  3.  It  is  interesting 

to  note  that  the  energy  solution  does  not  seem  to  be  converging  to  the 

natural  boundary  condition  at  x  =  0  and  a.  In  particular,  the  bending 

moment  M  is  increasing  as  it  approaches  the  boundary  rather  than  vanishing. 

Even  though  the  proper  bending  deflections  can  be  obtained  without  the 

natural  boundary  conditions  being  satisfied,  the  same  is  not  necessarily 

true  for  functions  such  as  bending  moments  which  depend  on  partial 

derivatives  of  the  deflection  function.  For  the  case  of  simple -supports  it 

is  very  difficult  to  chose  functions  in  conjunction  with  the  Ritz  method 

which  will  satisfy  the  natural  boundary  conditions.  In  particular,  due  to  the 

existence  of  D  ,  and  D  ,  cross -elasticity  terms,  the  moment  boundary 
16  Z6 

conditions  cannot  be  satisfied  by  a  set  of  assumed  functions  in  the  separable 


form  [  18] . 
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w  = 


(22) 


I  I  x«w 

m  n 


This  conclusion  is  based  on  the  assumption  that  all  desired  derivatives  of 

eq.  (22)  can  be  obtained  from  term-by-term  differentiation.  It  is  also 

anticipated  that  similar  difficulties  would  be  encountered  with  finite  element 

techniques  which  are  based  on  principles  of  minimum  potential  energy. 

The  convergence  of  eigenvalve  type  problems  is  qualitatively  similar 

to  the  maximum  deflection  results  as  illustrated  by  the  buckling  and  vibration 

o 

solutions  in  Tables  I  and  II  for  a  45  oriented  simply -supported  square  plate 
having  the  properties  of  eq.  (21)  .  Biaxial  compression  (K^  =  =  1,  =  0) 

and  pure  shear  (K^  =  =  0,  =  1)  are  both  shown. 


In  general,  for  a  laminated  plate  of  practical  construction  the  magnitude 

of  the  D  ,  and  D  .  bending  stiffness  terms  will  decrease  with  an  increase 
1  o  do 

in  the  number  of  plies.  Consider,  for  example,  a  symmetric  angle-ply 
laminate  (  +  8,  .  .  . ,  -  0,  -  0,  .  .  .,  +0  stacking  sequence)  .  The 
cross-elasticity  bending  stiffness  terms  are  of  the  form 
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h  Ql6(+0) 

4L 


D 


26 


h  Q26(+P) 
4L 


(23) 


while  the  remaining  bending  stiffness  terms  are  independem  of  the  number 
of  plies.  Thus,  the  effect  of  cros s -elasticity  bending  stiffness  terms  is 
less  severe  in  a  laminated  plate.  This  is  illustrated  in  Table  III  where  the 
deflection  and  bending  moment,  M  ,  at  the  center  of  a  4  layer  simply- 

X 

supported  ±  45°  (+45°,  -  45°,  -  45°,  +45°)  laminate  are  tabulated  for 
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increasing  values  of  M  and  N.  Unidirectional  ply  properties  are  those  of 
eq.  (21)  .  The  difference  between  the  Fourier  solution  and  the  Ritz 
solution  for  plate  deflection  is  much  less  than  in  the  case  of  a  singly  layer 
+  45°  plate.  It  should  be  noted,  however,  that  the  bending  moyfent  as 
determined  by  the  Ritz  method  displays  the  same  erratic  behavior  previously 
noted  by  Ashton  [10]  lor  a  homogeneous  anisotropic  plate.  In  particular, 
the  bending  moment  at  the  center  seems  to  be  os  'Ulatxug  about  a  value 
close  to  the  results  given  by  the  Fourier  solution  for  M  =  N  =  13. 
Furthermore,  the  amplitude  of  the  oscillation  s  .cms  to  become  larger  with 
increasing  values  of  M  and  N,  strongly  suggesting  that  M  does  not  converge. 

X 

Further  evidence  suggesting  such  a  conclusion  is  illustrated  in  Figure  5. 

The  distribution  of  M  across  the  center  of  the  plate  according  to  the  Ritz 

x 

solution  for  M  =  N  =  3  does  not  look  too  unreasonable  compared  to  the 
Fourier  solution.  For  M  =  N  =  7,  however,  the  Ritz  method  yields  results 
which  are  rather  horrifying. 

It  should  also  be  noted  that  the  Fourier  solution  yields  an  oscillating 
convergence  for  M  at  the  center  of  the  plate.  A  cursory  examination  of 

X 

Table  III,  however,  reveals  that  the  oscillations  damp  out  rapidly  with 
increasing  values  of  M  and  N. 

CONCLUSIONS 

Numerical  results  show  that  a  classical  Fourier  analysis  can  yield 
improved  solutions  for  homogeneous  and  symmetrically  laminated  anisotropic 
plates  having  strong  cross-elasticity  bending  stiffness  terms  compared  to 
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existing  solutions  obtained  by  energy  methods.  This  improvement  is  due 
to  the  fact  that  the  Fourier  solution  satisfies  both  the  geometric  and  natural 
boundary  conditions,  while  the  energy  method  in  conjunction  with  classical 
beam  mode  functions  satisfies  only  the  geometric  boundary  conditions. 
Convergence  of  the  Ritz  solutions  for  deflections,  buckling  loads,  and 
vibration  frequencies  seem  to  be  quite  slow  for  highly  anisotropic  plates 
in  which  all  four  edges  involve  natural  boundary  conditions.  It  is  quite 
questionable,  however,  whether  functions,  such  as  bending  moments,  which 
involve  derivatives  of  the  plate  deflection  converge  at  ail. 

The  Fourier  series  solution  appears  to  provide  a  conser/ative  solution 
for  the  results  presented.  In  particular,  bending  deflections  are  larger 
than  the  exact  solution,  while  buckling  loads  and  fundamental  vibration 
frequencies  are  low  compared  to  the  energy  solution  which  is  an  upper  bound. 
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TABLE  I 


Buckling  of  45 


o 


Simply -Supported  Plate, 


a  =  b 


2  3 

Nob  /Q -  Biaxial  Compression  (Kj  =  =  1 

i.  k3  =  0) 

1 - 

M 

N 

Fourier  Analysis 

Reference  7 

1 

1 

6.  763 

21.438 

3 

3 

8.  115 

13. 013 

5 

5 

8.  318 

11.565 

7 

7 

8.418 

11.0*0 

9 

9 

8.481 

11 

11 

8.521 

13 

13 

8.556 

N0b2/Q22h3  -  Uniform  Shear  (Kj  =  K2  =  0,  K3  =  1 ) 

1 

1 

12.121 

‘  JO 

3 

3 

13.309 

27.953 

5 

5 

,13.567 

18.022 

7 

7 

13.702 

17.  122 

9 

9 

13.786 

11 

11 

13.843 

13 

13 

13.884 
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TABLE  II 


Fundamental  Vibration  Frequency  of  45° 
Simply -Supported  Plate,  a  =  b 


o  b2  (p/Q!,2hV/2 

M 

N 

Fourier  Analysis 

Reference  7 

1 

1 

13.  277 

20.571 

3 

3 

14. 231 

17. 301 

5 

b 

14.415 

16. 737 

7 

7 

14.496 

16.428 

9 

9 

14.544 

11 

11 

14.577 

13 

13 

14.600 

j 

TABLE  III 


Bending  of  Simply -Supported  ±45°  Laminate 
(+45°,  -45°,  -45°.  +45°) ,  a  =  b,  q  =  qQ 


w(a/2,  b/2)  Q22h3/qob4x  103 

Mx(a/2,b/2)  /qQb2  x  102 

M 

N 

F  S 

R  S 

F  S 

R  S 

1 

1 

5.1838 

3.8310 

4. 7785 

4. 1854 

3 

3 

4.9676 

4.3850 

4.0025 

4. 2457 

5 

5 

4.8849 

4.3947 

4.2311 

3.6911 

7 

7 

4.8571 

4.4910 

4.1488 

4.2633 

9 

9 

4.8307 

4.1716 

11 

11 

4.8193 

4.1533 

13 

13 

4.8078 

4.1596 
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M 

FIGURE  2,  Maximum  Deflection  for 

I  ncreasing  Number  of  Terms 
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FIGURE  4.  Bending  Moment  M^x,  b/2),  Simply-Supported  Plate  Under  Uniform  Load, 
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x/a 

FIGURE  5.  Bending  Moment  M^x,  b/2),  Simply-Supported  Laminate 
(+45°,  -45°,  -45°,  +45°)  Under  Uniform  Lead,  a  »b 


